Computations of the shear stresses distribution experienced by passive particles as they circulate in turbulent flow: A case study for vWF protein molecules.

The stress distribution along the trajectories of passive particles released in turbulent flow were computed with the use of Lagrangian methods and direct numerical simulations. The flow fields selected were transitional Poiseuille-Couette flow situations found in ventricular assist devices and turbulent flows at conditions found in blood pumps. The passive particle properties were selected to represent molecules of the von Willebrand factor (vWF) protein. Damage to the vWF molecule can cause disease, most often related to hemostasis. The hydrodynamic shear stresses along the trajectories of the particles were calculated and the changes in the distribution of stresses were determined for proteins released in different locations in the flow field and as a function of exposure time. The stress distributions indicated that even when the average applied stress was within a safe operating regime, the proteins spent part of their trajectories in flow areas of damaging stress. Further examination showed that the history of the distribution of stresses applied on the vWF molecules, rather than the average, should be used to evaluate hydrodynamically-induced damage.

Distribution and history of extensional stresses on vWF surrogate molecules in turbulent flow.

The configuration of proteins is critical for their biochemical behavior. Mechanical stresses that act on them can affect their behavior leading to the development of decease. The von Willebrand factor (vWF) protein circulating with the blood loses its efficacy when it undergoes non-physiological hemodynamic stresses. While often overlooked, extensional stresses can affect the structure of vWF at much lower stress levels than shear stresses. The statistical distribution of extensional stress as it applies on models of the vWF molecule within turbulent flow was examined here. The stress on the molecules of the protein was calculated with computations that utilized a Lagrangian approach for the determination of the molecule trajectories in the flow filed. The history of the stresses on the proteins was also calculated. Two different flow fields were considered as models of typical flows in cardiovascular mechanical devises, one was a Poiseuille flow and the other was a Poiseuille-Couette flow field. The data showed that the distribution of stresses is important for the design of blood flow devices because the average stress can be below the critical value for protein damage, but tails of the distribution can be outside the critical stress regime.

GPU-based detection of protein cavities using Gaussian surfaces.

BACKGROUND: Protein cavities play a key role in biomolecular recognition and function, particularly in protein-ligand interactions, as usual in drug discovery and design. Grid-based cavity detection methods aim at finding cavities as aggregates of grid nodes outside the molecule, under the condition that such cavities are bracketed by nodes on the molecule surface along a set of directions (not necessarily aligned with coordinate axes). Therefore, these methods are sensitive to scanning directions, a problem that we call cavity ground-and-walls ambiguity, i.e., they depend on the position and orientation of the protein in the discretized domain. Also, it is hard to distinguish grid nodes belonging to protein cavities amongst all those outside the protein, a problem that we call cavity ceiling ambiguity. RESULTS: We solve those two ambiguity problems using two implicit isosurfaces of the protein, the protein surface itself (called inner isosurface) that excludes all its interior nodes from any cavity, and the outer isosurface that excludes most of its exterior nodes from any cavity. Summing up, the cavities are formed from nodes located between these two isosurfaces. It is worth noting that these two surfaces do not need to be evaluated (i.e., sampled), triangulated, and rendered on the screen to find the cavities in between; their defining analytic functions are enough to determine which grid nodes are in the empty space between them. CONCLUSION: This article introduces a novel geometric algorithm to detect cavities on the protein surface that takes advantage of the real analytic functions describing two Gaussian surfaces of a given protein.

Turbulent plane Poiseuille-Couette flow as a model for fluid slip over superhydrophobic surfaces.

In this study, plane Poiseuille-Couette flow is simulated as a model for specified streamwise slip on one of the channel walls. The relative velocity between the two walls is set to be 1, 2, and 4 in viscous wall units. This is equivalent to the presence of a superhydrophobic surface at one of the channel walls that causes fluid to slip on the boundary. The results show that the streamwise slip forces turbulence in the near-wall region to tend towards a limiting one-component state. This leads to the suppression of small scale turbulence and laminarization close to the wall and then to drag reduction. The selective weakening of the streamwise vorticity close the wall and the observed decrease of turbulence kinetic energy production can then be considered as a consequence of this effect. Changes in the coherent structures, including a decrease of sweep events and increase of ejection events close to the wall where slip occurs, are also observed.